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VII. — On the Theory of the Perturbations of the Planets. By James Ivory, A.M. 
F.R.S. Instit. Reg. Sc. Paris. Corresp. et Reg. Sc. Gottin. Corresp. 

Read January 19, 1832. 

1 HE perturbations of the planets is the subject of reiterated researches by all 
the great geometers who have raised up Physical Astronomy to its present 
elevation. They have been successful in determining the variations which the 
elements of the orbit of a disturbed planet undergo ; and in expressing these 
variations analytically, in the manner best adapted for computation. But the 
inquirer who turns his attention to this branch of study will find that it is 
made to depend upon a theory in mechanics, which is one of considerable 
analytical intricacy, known by the name of the Variation of the Arbitrary 
Constants. Considerations similar to those employed in this theory were 
found necessary in Physical Astronomy from its origin; but the genius of 
Lagrange imagined and completed the analytical processes of general appli- 
cation. In a dynamical problem which is capable of an exact solution, such 
as a planet revolving by the central attraction of the sun, the formulas con- 
structed by Lagrange enable us to ascertain the alterations that will be in- 
duced on the original motions of the body, if we suppose it urged by new 
and very small forces, such as the irregular attractions of the other bodies of 
the planetary system. General views of this nature are very valuable, and 
contribute greatly to the advancement of science. But their application is 
sometimes attended with inconvenience. In particular cases, the general 
structure of the formulas may require a long train of calculation, in order 
to extricate the values of the quantities sought. It may be necessary for at- 
taining this end to pass through many differential equations, and to submit to 
much subordinate calculation. The remedy for this inconvenience seems to 
lie in separating the general principles from the analytical processes by which 
they are carried into effect. In some important problems, a great advantage, 

2 c 2 
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both in brevity and clearness, will be obtained by adapting the investigation 
to the particular circumstance of the case, and attending solely to the princi- 
ples of the method in deducing the solution. It may therefore become a ques- 
tion whether it be not possible to simplify physical astronomy by calling in 
the aid only of the usual principles of dynamics, and by setting aside every 
formula or equation not absolutely necessary for arriving at the final results. 
The utility of such an attempt, if successful, can hardly be doubted. By ren- 
dering more accessible a subject of great interest and importance, the study 
of English mathematicians may be recalled to a theory which, although it 
originated in England, has not received the attention it deserves, and which 
it has met with in foreign countries. 

The paper which I have the honour to submit to the Royal Society, contains 
a complete determination of the variable elements of the elliptic orbit of a dis- 
turbed planet, deduced from three differential equations that follow readily 
from the mechanical conditions of the problem. In applying these equations, 
the procedure is the same whether a planet is urged by the sole action of the 
central force of the sun, or is besides disturbed by the attraction of other 
bodies revolving about that luminary ; the only difference being that, in the 
first case, the elements of the orbit are all constant, whereas in the other case 
they are all variable. The success of the method here followed is derived from 
a new differential equation between the time and the area described by the 
planet in its momentary plane, which greatly shortens the investigation by 
making it unnecessary to consider the projection of the orbit. But the solution 
in this paper, although no reference is made to the analytical formulas of the 
theory of the variation of the arbitrary constants, is no less an application of 
that method, and an example of its utility and of the necessity of employing it 
in very complicated problems. 

1. If S represent the sun and P, P two planets circulating round that lumi- 
nary, it is proposed to investigate the effect of the attraction of P to disturb 
the motion of P and to change the elements of its orbit. We here confine our 
attention to one disturbing planet ; for there is no difficulty in extending to 
any number, the conclusions that shall be established in the case of one. 

The positions of the planets P and P may be ascertained as usual by the 
rectangular coordinates x, y, z and x', y', z' ; x, y, x' } y being contained in a 
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plane passing through the origin of the coordinates placed in the sun's centre ; 
and z, z' being perpendicular to the same plane. 

Further, let M, m, m! denote the respective masses of S, P, P' ; r and r' the 
distance of P and P' from S, and § the distance between the two planets ; then, 

putting {ij = M + m, the direct attraction between S and P will be -% -, and 

the resolved parts of this force, acting in the respective directions of x, y, z, and 
tending to diminish these lines, will be 

(J.X ft,1/ fi.Z 
~^W) ~^3) ~^3' 

The planet F attracts S with a force = -p%, of which the resolved parts are, 

m' or? m'y' wf z' 



jJS } r '3 ) y/3 • 

The same planet P' attracts P with the force t, of which the partial forces 

m! {x' — x) m' (y — y) w! (z' — z) 

Were S and P attracted by F in like directions with equal intensity, the rela- 
tive situation of the two bodies would not be changed, and the action of P 
might be neglected : but the attractions parallel to the coordinates being un- 
equal, the differences of these attractions, viz. 

m' {x f — x) m'x 1 m'(y'—y) m'y' m' (z' — z) m' z 1 
3 r t$ ) 15 ~ r rs 5 3 r '3 3 

are exerted in altering the place of P relatively to S. These last forces increase 
the coordinates x, y, z ; and, therefore, they must be subtracted from the 
former forces which have opposite directions, in order to obtain the total 
forces acting in the directions of the coordinates and affecting the motion of 
P relatively to S, viz. 

m x m! (x' — x) m' x 



7.3 „s 


' 1" f J3 : 


my m' (y' — y) 


m'y' 


mz m! (z' — z) 


m' z' 
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But, if dt represent the element of the time supposed to flow uniformly, the 

actual velocities with which the coordinates increase are, Ji> Tp di' an( * ^ e 

increments of these velocities, -j£ t -j^, -j^r, are the effects produced by all 

the forces that urge the planet. Equating now the forces really in action to 
the measure of the effects they produce, and observing that the two equivalent 
quantities have been estimated in opposite directions, we obtain the following 
equations for determining the place of P relatively to S at any proposed instant 

of time, 

ddx ix./x m' (x' — x) m' x' 



~r ~s — 



df "•" r 3 ~ g 3 " r' 3 > 

ddy , f*J/ _ m'jy 1 -y) _ ^£ 
d t* ~r r 3 ~~ f r 13 ' 

ddz pz _ m' (z f — z) __ irH 
df + ~r» ~~ g 3 r' 3 ' 

If we now assume 



„ __ to[ f i __ xd +y.y' + zz ' 1 

**' x I • (*' - *)■ + (y -yf + (*'^p (^ + y a + z'*f S ' 



, , , . , ,.„ A . , dR dR dR 

it will be found that the partial differentials, \b X -j~ x , y> X jr, [* X -j^, are 



respectively equal to the quantities on the right sides of the last equations, 
that is, to the disturbing forces tending to increase the coordinates w, y, z. 
These equations may therefore be thus written, 
ddx x d R 

i*y. j_ a. ~ 15 v (A) 

p.dP + r 3 — dy> }■••■■■■•■■■■ \*) 

ddz z dR 

JTd? + t* ~ ~dz> 

If it be asked, What notion must be affixed to the symbol pelt 2 }, it will be 

recollected that 4- is the attraction between S and P at the distance r ; and if 

we suppose that P describes a circle, of which unit is the radius, round S, the 

centripetal force in the circle will be |r or p ; and the velocity with which P 
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moves in the circle will be proportional to ^/p. Thus the algebraic quantities 
ttjp' and dt,J^ represent the arcs of this circular orbit, which are described 
in the times t and dt. 

It is requisite in what follows to transform the coordinates x, y, % into other 
variable quantities better adapted for use in astronomy. Let X and k' denote 
the longitudes of the planets P and F reckoned in the fixt plane of x y, and 
s and s the tangents of their latitudes, that is, of the angles which the radii 
vectores r and r' make with the same plane : then, 

r cos A , _ r 1 cos A' 



a/1+s*' ~ * V'l + s' 2 ' 
r sin A f _ r' sin a' 

y ~ vT+T 3 ' y ~ a/I + s' 2 ' 

r s , r' s 1 



a/1+s 3 ' a/1 + s' 2 ' 

In the transformations alluded to, the quantities -j^, -j-, -j- must be ex- 
pressed in the partial differentials of R relatively to the new variables r, X, s ; 
and it will conduce to clearness of method if these calculations be dispatched 
here. We have the equation, 

dR __ dR dr dR d\ dR ds^ 
dx dr ' dx ' d/\ ' dx ' ds ' dx ' 

and having computed the differentials jj, j^,, j^ from the formulas 



r = Jx* + y* + z\ tan.X = f, s = 



>/x* + f> 

the substitution of the results will make known the expression of -j- . By the 
like procedure the values of -j- and -j^ will be found 



dR _dR 
dx ~ dr ' 


cos A d R 
a/1+s* d\ 

sin A dR 

*/\ + s * + dx ' 

s 


sin A a/1 + s s dR 
r ~ ds ' 


cos Asv'l + s a 
r 


dR dR 

dy dr ' 


cosAy'l + s 3 dR 
r ds 

dR 

+ Is 


sin As a/1 + s 2 
r 


dR dR 


V\ +s 3 


dz dr ' 


a/I +S 3 


r 



>■ • • (B) 
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The new partial differentials of R represent the disturbing forces reduced to 
new directions. By combining the formulas (B), we get 

d R d R cos \ ^ R sin A dR s 

~dr z ~~dx' v'l +s 9 "*" ^y " V'l + ** ^ ' V'l + s 2 ' 

and it will readily appear that the coefficients of -j^, -j~, -j^ are the re- 
spective cosines of the angles which the directions of the forces make with r -, 
so that -j- is the sum of the three partial forces that urge the planet from the 



sun. In like manner it may be proved that -^- . r + s is the disturbing 

force perpendicular to the plane passing through the sun and the coordinate z, 

that is, to the circle of latitude ; and that -jj . — — is the force acting in the 

same plane perpendicular to r, and tending to increase the latitude. 

2. If the equations (A), after being multiplied by 2 dx, 2dy, 2dz, be 
added together, and then integrated, we shall get this well-known result, 

da* + df + dz* 2 l _ 9/v ,t, m 

in which — is the arbitrary constant, and the symbol d' R is put for 

dR , , dll . , dR , 
-^dx + ^dy+^dz; 

that is, for the differential of R, on the supposition that x, y, z, the coordi- 
nates of the disturbed planet, are alone variable. If we conceive that R is 
transformed into a function of the other quantities r, X, s, we shall therefore 

have 

7 , „ dR j , dR j dR , 
d'R= in :dr + in: dX+ in ds. 

Supposing that the radius vector r, at the end of the small interval of time 
dt, becomes equal to r + dr, and that dv expresses the small angle contained 
between r and r -{- dr,we shall have 

d r i _j_ r i dip. = d x 2 _J- dy2 -|_ d z i ; 

for each of these quantities is equal to the square of the small portion of its 
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orbit which the planet describes in the time d t. The last equation may there- 
fore be thus written, 

J^ + ^_A + _L = 2 /v R (2) 

fKlr ' par r ' a J y ' 

The double of the small area contained between the radii r and r + dr, is 
equal to r^dv, and as x, y, z and x + dx, y + dy, z + dz, are the coordi- 
nates of the extremities of the radii, the projections of the area upon the planes 
of xy, xz, yz, are respectively equal to 

xdy — ydx, xdz — z dx, ydz — zdy: 
wherefore, according to a well known property, we shall have, 

r*dv 2 (xdy — ydxf (xdz — zdxf (yds; —zdy) 2 

yd¥ ~ jTIF •" pdf "l pdr : 

and the differential of this equation, dt being constant, may be thus written, 

7 r*dv* „ . „ , „ , „, /dxddx-\-dyddy+dssddz\ 

d --^ = 2 ^ 2 +y 2 + z ) • \ — — Ur ) 

-2(xdx+ydy + zdz) . ( ****+** '2 +1***). 

Now, substitute the values of the second differentials taken from the equations 
(A), and we shall obtain, first, 



dx ddx + dyddy + dssddz djl , , dR , d_R , djr ^ 

'~£d? — dx ax + dy d y '+ dz dZ ~ r* — d K 

and, secondly, 

xddx+yddy + zddz dR. <2 R 

^dr ~~~ ~~dx X ~T~ ~dyy "•" 



dr 



— _ J- — -~ __ _L 
dz r dr r 



wherefore, since x 2 + y 2 + ^ 2 = r 2 and xdx-{-ydy-\"zdz = rdr, the fore- 
going differential equation will become by substitution, 



or, which is equivalent, 



7 r*dv* „ „ / 7 ,„ rfR , \ 



'.££ = "(£" + £«)• 
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By integrating, 

r 2 dv = hdt,J[h, 



h* = h * + 2fr*(d'R-?£dr), 



dR 

~j — i 
ar 

dR ,. . dR 



> (3) 



the constant h being equal to — - when tf = 0. 

Further, let the first of the equations (A) multiplied by y be subtracted from 
the second multiplied by on ; then 

d.(xdy-ydx) dR dR 

pde dy x d x y • 

and, by converting the quantities in this equation into functions of r, X, s, 

7 / r* d\ \ 

a - \i +5 3 ' dty-y,) _dR, 

and by multiplying both sides by 2 . - ■ - ^ . rfx, 

, / » ,a <?*■ \2 _ r 3 dR 



and, by integrating, 

r % d\ 

1 +s 2 



= h' dt^/p, 



1 



J> (4) 

^ = V* + 2/tt?;^J 

# ' being a constant. 

The equations that have been investigated, which are only three, the first 
and second being one equation in two different forms, are sufficient for deter- 
mining the place of a planet at any proposed instant of time, whether it revolves 
solely by the central force of the sun, or is disturbed by the irregular attrac- 
tions of the other bodies of the system. The second and third equations 
ascertain the form and magnitude of the orbit in its proper plane, and the 
place of the planet ; the fourth equation enables us to find the angle in which 
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the plane of the orbit is inclined to the immoveable plane of x y, and the posi- 
tion of the line in which the two planes intersect one another. 

3. We begin with the more simple case of the problem, when the planet is 
urged solely by the central force of the sun. On this supposition, there being 
no disturbing forces, we must make R = in the equations of the last §. 
By the formulas (3) and (4), we have, 

r 2 d v — hdt */(*>, 
Yj—z.dX^h'dttJp; 

and h, k', are constant quantities. Now— r====== is the projection of r upon 

r 3 
the plane of x y ; and the area j-t~^ >dX is the projection of the area r 2 dv 

upon the same plane ; wherefore, if i denote the angle of inclination which the 
plane containing the radii vectores r and r-\-dr, has to the plane of x y, we 
shall have 

C0Sl== r'dv =T* 

which proves that a plane passing through the sun's centre and any two places 
of the planet infinitely near one another, has constantly the same inclination 
to the immoveable plane of xy. And it further proves that the planet moves 
in one invariable plane ; for, unless this were the case, the areas described 
round the sun in any consecutive small portions of time, could not constantly 
have the same proportion to their projections upon the plane of xy. 

The orbit in its proper plane will be determined by the equations (2) and 
(3), viz. 

dr* . i*d<d 



+ n.rli? r + n ~ °5 



p.df~ pd? r ~ a 

r 2 dv = hdt *Jp> 

a and h being arbitrary quantities. By exterminating dt„Jp from the first 
equation, 

2 d2 
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7.2 
by multiplying all the terms by — , and adding 1 to both sides, 

} JL d7& .1/1 r \ 2 — h * 

a ' r*dv 2 ~T~\ l a) ~ l a ' 
and by introducing the new quantity e 2 , 

2=1-- 

a 
This equation is solved by assuming 

dr e sin 3 

r d v 1 + e cos fl 



r . , cos + e 

a 



~ " e X 1 + ecosfl' 



the arc 6 remaining indeterminate. For, if the assumed quantities be sub- 
stituted, the equation will be verified, and the arc 6 will be eliminated. In 
order to determine 6, let the second of the formulas be differentiated, and equate 

— to the like value in the first formula ; then, 
r 

dv = dd', and v — w = 6. 
The nature of the orbit is therefore determined by these two equations, 

dr e sin (v — is) 

r dv 1 + e cos (» — ■or)' 

a(l-<?) 



r = 



1 + e cos (v — ts) ' 



the first of which shows that the two conditions -77 = 0, and sin (v — nr) = 0. 

must take place at the same time ; so that to- is the place of the planet when 
its distance from the sun is a minimum = a (1 — e), or a maximum = a (I + e) : 
and the second proves that the orbit of the planet is an ellipse having the sun 
in one focus ; a being the mean distance ; e the eccentricity ; and v — w the 
true anomaly, that is, the angular distance from the perihelion or aphelion ; 
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from the perihelion if e be positive, and from the aphelion if the same quan- 
tity be negative. 
It must however be observed that the preceding determination rests entirely 

on the assumption that, in the equation e 2 = 1 — — , the quantity — is posi- 
tive and less than unit. Without entering upon any detail, which our present 
purpose does not require, all the possible cases of the problem will be suc- 
cinctly distinguished by writing the equation in this form, 

h* 



1 +e 



= (1 — e) X a. 



The quantity on the left side being essentially positive, the two factors on the 
other side must both have the same sign. If they are positive, the orbit will 
be an ellipse ; if they are negative, and consequently e greater than unit, the 
curve described by the body will be a hyperbola ; and it will be a parabola, 
when e = 1, and a and 1 — e pass from being positive to be negative, at which 
limit the equation will assume this form, 

— = X <*3. 

¥■ 
In all the cases Y~Z~ e * s tne perihelion distance. 

The nature of the orbit being found, we have next to determine the relation 
between the time and the angular motion of the planet. For this purpose we 
have the equation, r 2 dv = hdt s /(^, from which, by substituting the values 
of r and h, we deduce 

#4- M + e cos (v — ■en-) ) 

Let — -t = n ; then, by integrating, 



nt -\- s — ts — 



(l -fpjlv 



J (\ + e cos (v — ■&) ) 



the quantity under the sign of integration being taken so as to vanish when 
v — zs = 0, and s being a constant quantity. The mean motion of the planet 
reckoned from a given epoch, is equal to nt + s; and the mean anomaly, to 
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n t + s — w, the true anomaly being v — -ar. The equation may be put in this 
form, 

J 1 + e cos {v — w) 1 + e cos (v — w) 

and, if we assume 

// 1 — e 3 . sin (v ■*- «r) cos (p — ot) + g 

sina ~ 1 + ecos(w - v) > C0S M _ 1 + ecos(w-w) ; 

we shall find, 



/ * Vi —f.dv 
J 1 + e cos (» — ■ « 



+ e cos (p — «■) 
so that we readily arrive at these results, 
nt -\- & — is = u — e sinu, 

r = l +<? cos(t,-«r) = "(! ~ ^OS«), 



w — -us u f\ a. e 

tan -y- = tan ... " -./ \ T * 



W^i 



These last are the formulas that occur in the solution of Kepler's problem, 
the arc u being the anomaly of the eccentric. Having found the expression of 
the eccentric anomaly in terms of the mean anomaly from the first of the for- 
mulas, we thence deduce the true anomaly v — w, and the radius vector r, for 
any proposed instant of time. The analytical solution of these questions is 
omitted ; the sole intention of treating here of the motion of a planet circu- 
lating by the central force of the sun, being to elucidate the investigations that 
are to follow respecting the orbit of a disturbed planet. 

The purposes of astronomy require further that the motion of the planet in 
its orbit be connected with the longitudes and latitudes estimated with regard 
to the immovable plane of xy. The orbit being supposed to intersect the im- 
movable plane, and the angle of inclination being represented by i, let N stand 
for the longitude of the ascending node, and P for the place of the same node 
in the plane of the orbit and reckoned from the same origin with the true 
motion v : then v — P, or the distance of the planet from the node in the 
plane of the orbit, is the hypothenuse of a right-angled spherical triangle, one 
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side of which is the arc X — N in the immovable plane, and the remaining side 
is the latitude having s for its tangent : wherefore we have 

tan (X — N) = tan {v — P) cos i, 

s = tan i sin (X — N). 

The first of these equations enables us to compute X when v is given, and con- 
versely ; by means of the second, the latitude is found. The practical calcula- 
tions are much facilitated by expressing the quantities sought in converging 
serieses : but the discussion of these points is beside our present purpose. 

4. We now proceed to investigate the effect of the disturbing force of the 
planet P in altering the orbit of P. For this purpose we have the equations 
(3) and (4), viz. 

r 2 dv = hdt *Jiij, 

dX = h! dt dp; 



l +s 2 



of which the first is the expression of the small area described round the sun 
by the planet in the time d t, and the other is the projection of that area upon 
the immovable plane of xy. Wherefore, if i denote the angle of inclination 
which the plane passing through the sun and the radii vectores r and r + dr, 
has to the plane of xy, we shall have 

-.2 



C0SU— r*dv ~~ h ' 

and, as h' and h vary incessantly by the action of the disturbing forces, it 
follows that the momentary plane in which the planet moves is continually 
changing its inclination to the fixed plane. Let i' be the value of i when t = ; 

h' 

then cos i' = -f ; and, by the formulas (3) and (4), we shall have, 

fin 






o ™ a " ~ry ' ■ dx • l +s 3 ' 
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and hence, in consequence of what has been shown, 

„. A» . -. A" 3 F 2 

cosn = ^5- ; sin 2 i = -^ ; tann = -p. 

Let the momentary plane of the planet's orbit, that is, the plane passing 
through r and r + dr, intersect the immovable plane of x y, and put N for the 
place of the ascending node: then s and s + ds will be the tangents of the 
latitudes at the distances X — N, and X + dx— N from the node : and, i being 
the angle contained between the two planes, we shall have, 
s = tan i sin (X — N), "") 

ds . . , XTX r • ( 5 ) 

■j^ = tan i cos (X — N). I 

By adding the squares of these equations, 

ds* . h" 3 

* ~*~ dl? = ^ an l == Tfi ' 

by differentiating, making J X constant, 

h»dk"~tidh'(s* + ~\ 



Hi£+*)= 



ds 

and, by substituting the values of h" dh" and h' d h', 

dds t^ CdR I dR ds") 

tf A* + * — h' 3 • \ds ~ I + s 2 • dx • dx]- 

Since i is variable in the equations (5), it is obvious that N, or the place of 
the node, must likewise vary. By combining each of the two equations with 
the differential of the other, these results will be obtained, 

d . tan i . _ . dN^ . ._ T . 

= dx sin (X — N) — -^- tan ^ cos (X — N) 

dds d.tani XT . , tfN . . . XTX 

dtf + * = ~3X* C0S ^ _ N ) + U tan * Sm ( X ~ N )' 

from which we deduce, 

J i = cos 2 i cos (X — N) . < -ttj + s > . J X, 

, , T cos i sin (A — N) ( dds , 1 , . 
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and, by substituting the value of -j^ + s, 

and, observing that r 2 d X = (1 + * 2 ) h' d t ^/^ = (1 + * 2 ) k cos idt ,Jp, we 
finally get, 

« = «Je«(X-M).{(l + ^£-g.£}.«{tf] 

*N=^>.{(.+^-^.£}.^ j'' (6) 

These equations determine the motion of the node in longitude, and the varia 
tion in the inclination of the orbit. They are rigorously exact, and may be 
transformed in various ways, as it may suit the purpose of the inquirer. 

We proceed now to investigate the motion of the planet round the sun. 
For this purpose we have the equations (2) and (3), viz. 

pdt 2 ^ pd* r T a — *J « «■> 

r 2 dv = hdt /s/p. 

And first, as the small arc dv contained between the two radii r and r + dr, 
continually passes from one plane to another, it is requisite to inquire what 
notion we must affix to the sum v. The momentary plane of the planet's 
motion, in shifting its place, turns upon a radius vector ; and if we suppose a 
circle concentric with the sun to be described in it, and to remain firmly at- 
tached to it, the differentials dv will evidently accumulate upon the circum- 
ference of the circle, and will form a continuous sum, in the same manner as 
if the plane remained motionless in one position. The arc v is therefore the 
angular motion of the planet round the sun in the moveable plane, and is 
reckoned upon the circumference of the circle from an arbitrary origin. 

In the first of the foregoing equations a is an arbitrary constant, and I shall 
put, 



i=i-.yvH, 



so that we shall have 

dr* r*dv* _2 J. 

ltdt* + pdf ~ r + a ~"° 5 

r 2 dv = hdtfijp; 
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which are different from the corresponding equations in the last section in no 
respect, except that here h and a are both variable, whereas in the other case 
they were both constant. Treating these equations exactly as before, we first 
get by exterminating dt/^/JL, 

jn dr* J? z_ j_ 

r % 
then, by multiplying all the terms by -^ and adding 1 to both sides, 



a ' t»d* + \ 1 "* a) — 1 ~- a ' 



from which we deduce 

e 2 = i _ A 2 

a ■ 



^^S . j, v 2 



(»-Oi?fc+'0-f)=*- 

The last equation is solved by the same assumptions as before, viz. 



dr e sin 



rdv 1 + e cos 0' 

r cos + e 

1 *_ — p v ' — • 

a -" e * 1 + ecos0' 

but it must be recollected that in these formulas, a and e are both variable. 
By differentiating the expression of r, viz. 

— a (1 - e a ) y 

r ~~ 1 + e cos 9 ~~ 1 + e cos 0' 

we get, 

<Zr esinfi.rfS £<2A cos0. <?£ 

r 1 + e cos ' h ~~ I + e cos ' 

d t 
and by equating this expression to the value of— taken from the first formula, 

and reducing, we obtain, 

2 dh 
e (dv — d&) sin 6 + cos 0.de= ~j- (1 + p cos 6). 

It appears therefore that v — 6, or ts, is a variable quantity; and the formulas 
that determine the elliptic orbit, and the variation of w, are as follows : 
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d r e sin (v — «r) 



r d v 1 + e cos (v — •sr) 

— a(l -e*) ___ g 

""" 1 + e cos (» — w) 1 + e cos (v — tsf 

Qdh ( \ 

ed7ffsin(v — 7s)+de.cos(v — 7s)==-j-{l-{-ecos(;v—&)j. ... (7) 

It is obvious that this last formula is tantamount to the equating to zero of 
the differential of r relatively to the variables, h, e, &, or a, e, & ; it may there- 
fore be thus written, 

T-rfa+ j-de-\- t— dn = (8) 

da ' de ' aw v ' 

The equations that have been investigated, enable us to deduce from the 
disturbing forces the variable elements of the ellipse that coincides momen- 
tarily with the real path of the planet ; a being the mean distance, e the eccen- 
tricity ; m the place of the perihelion, and h 2 the semi-parameter. We have 
next to find the relation between the time and the angular motion in the 
variable orbit. This will be accomplished by means of the equation r 2 dv 
= hdt J^ ; from which we obtain, by substituting the values of r and h, 

dt V> (1 -e^.dv 

a* (l + ecos(w — to)) 

The integral C 3 ^ , supposed to commence with the time, is the mean 

motion of the planet : when there is no disturbing force, a being constant, the 
mean motion is proportional to the time and equal to n X t ; but the action of 
the disturbing forces, by making a variable, alters the case, and requires the 
introducing of a new symbol £ to represent the mean motion. Thus we have 

t— C~ X dt; dl = ~ -ci- 

* J & + ^cos(»-w)) 

The value of t, cannot be obtained directly by integration on account of the 
variability of e and vs. Let/(w — -m, e) express that function of the true ano- 
maly which is equal to the mean anomaly in the undisturbed orbit ; that is, 
suppose, 

2 E 2 
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J V J (l +ecos(y-OT)) 

the integral being taken on the conditions that it vanishes when v — sr = 0, 
and that <? and w are constant. If now we make e and w variable, we shall 
have, 

d.flv — ^e) , , d.f(v — w,e) , , d,f(v — w,e) 7 7 /./ \ 

dv dv +~ de de + dv - d* = a.f(v-*,e)- 

But the partial differential relatively to v, is no other than the expression of 
<Z£: wherefore, 

d£ + JK de de+ JK d ^ 'd* = d.f(v-*,e). 

By introducing a new symbol this equation may be separated into the two 
which follow, 

dZ,-^- ds — dtf = d .f(v — &, e), 

, d . f(t) — ix.e), , d. f(v — ■ar.e) , 

ds-d*=: J( de U e+ J \^ rfg. 

In the integral 

£-j-S_7tf — f{v—73,e), 

g + g is the mean motion of the planet reckoned from a given epoch, e however 
representing a quantity that varies incessantly by the action of the disturbing 
forces, the amount of the variation being determined by the second formula 
in which the value of s alone has not been previously ascertained. The mean 
anomaly of the planet is £ + s — ■& ; and the integral shows that there is the 
same finite equation between the mean and the true anomalies in the disturbed 
orbit as when there is no disturbing force. It follows therefore that, in both 
the cases, the true anomaly, the true motion of the planet, and the radius 
vector, are deducible from the mean anomaly by the same rules and by the 
solution of Kepler's problem. 

In order to find the value of the new variable i, it is necessary to eliminate 
the differential coefficients from its expression. Differentiating relatively to e 
and sr, we shall get, 



d.fjv — ot, e) 
de 



, /-.21 cos (v— si-) + 3e + e" cos (v—&) , 

= — ^l\ — e 2 . / • 7ST • dv, 

v L e -J (l + ecos(w-w)) 
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d .f(v — 1st, e) _ , „!- /» 2 e sin (v — w) . d v 



and, by integrating, 



<f .f(v — ot, e) sin (g-w)^!-^ sin (v — ■a) >/ 1 — e* 

Te 1 + ecos (v - w) (i + ecos („ _ w )y " 



( 2 +■ e cos (w — otm sin (v — is) »/ \ — e 
(i + ecos (v — wn 



d.f(v-«r,e) _ (l-<?)* 

These values being substituted in the foregoing formula, we shall find this re- 
sult, after dividing all the terms by the coefficient of dw, 

(l + e cos (v — w)Y (2 + e cos (v — &)) sin (w — cr) 

- — . (ds — dvx) = — - : -r ««- da, 

or, more concisely, 

02 1/1 3 (^2 + £?cos(» — -Br) )sin(w— «r) 

a ^g .(rf«-rf*)=»-i-^ \_J '.de-d*.. . .(9) 

From the equation between the mean and the true anomalies we deduce, 

<1> representing a function of the mean anomaly Z, + £ — OT "• an d as the differ- 
entials of £ and v are independent of the differentials of s, e, and ™, we shall 

have, 

dw , . dv , , dv , . ,, _v 

_ rfg + _J e + _^ = (10) 

Now, 

and, because, O is a function of s — *?, 

d.$ d.# ., dv . dv . 

' 1f7 = - SZ : consequently, ^ + ^ = 1 . 

The equation may therefore be thus written, 



214 MR. IVORY ON THE THEORY 

But, v being a function £ + *> it follows that, 



dv__ dv _ s? V 1 — e* , 
rfs ~~ if? r 2 ; 

and thus it appears that the equation we are considering is identical with the 
formula (9) : from which we learn that, 

d v (% + e cos (v — otm sin (v — w) 

__ _ ___ . 

It remains now to say a word about the longitudes and latitudes of the 
planet reckoned on the immoveable plane of xy. The variable quantities N 
and i denote the longitude of the ascending node, and the inclination of the 
orbit, in respect to the fixed plane : let P represent the place of the same node 
on the moveable plane of the planet, this arc being reckoned from the same 
origin as the true motion v : then, because the momentary plane in which the 
planet moves, in taking a new position, turns about a radius vector, it is ob- 
vious that, if d N be the motion of the node in the fixed plane of xy, cos i X dN 
will be its motion in the variable plane of the orbit. Wherefore we have, 

dV = cos i x ^N, and P —fcos i . dN, 

a constant being supposed to accompany the integral. This being observed, 
it is obvious that the same equations as in the case of the undisturbed orbit, 
will obtain between the quantities under consideration, viz. 

tan (X — N) = cos i tan (v — P), 

s = tan i sin (X — N). 

The foregoing investigations prove that the motion of a disturbed planet 
may be accurately represented by a variable ellipse coinciding momentarily 
with the real path of the planet. The variations, in the magnitude, the form, 
and the position of the ellipse, have been expressed by equations that depend 
upon the disturbing forces. A new inquiry presents itself: to exhibit the 
differentials of the elements of the variable orbit in the forms best adapted 
for use. 
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5. The expressions of the coordinates a?, y, z, in terms of the variables 
r, X, s, are as follows : 

r cos X r sin X _ r s 



V' 1 + s 2 ' y V 1+ s 3 ' Vl+s 8 " 

and, if we write X — N + N for X, we shall get, 

f cos (X — N) , T sin (X — N) . ._, ") 

^ = r 4~vTTF cosN --^TTT smN j' 

{sin (X — N) XT . cos (X — N) . - , "> 

But t; — P in the plane of the planet's motion is the hypothenuse of a right- 
angled spherical triangle of which X — N is one side, s the tangent of the 
other side, and i the angle opposite to this latter side ; and from these consi- 
derations we get 

cos (X — N) , „, sin (X — N) . . „. . , 

— ^==^ = cos(t;-P), J i+s J = sm(v-P)cosi, and 

sin (v — P) sin i : 



a/i + < 
wherefore we have these values of the coordinates, 

x = r . {cos (w — P) cosN — sin (v — P) sin N cos i} 
y = r. {sin (v — P) cos N cos i + cos (v — P) sin N} 
z = r . sin (v — P) sin i. 

The radius vector r is a function of v, a, e, sr, viz. 

— a(l-e a ) 

1 + e cos (w — ■or) " 

and thus the coordinates x, y, z, are functions of v and the five elements a, e, 
tg, N, z ; for P is no independent quantity, since it varies with N. In order to 
abridge we may write X, Y, Z for the multipliers of r in the foregoing expres- 
sions of x, y, z ; so that 

x = r X X, y = r X Y, z = r X Z. 

Now, on account of the equation (8) we have 
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dx 

d a 



, , dx 7 , dx 7 (dr , , dr 7 , dr 7 1 _ r _. 

da + Te de + dZ d W = \T & da + Te de + dZ d ") X X==0 



and, in like manner, 



-r- dz,-\--r de + -r~ dsr = 0, 

aa ' a« ■ aw ' 



rfs 



ds 



l- a d ^ + JZ de + 



dz 



da. 



de 



dx 



0. 



Further, we have, 



* -7xt , dx ,. C /dX . . dX\ dX,.\ 

x dN + dJ dz = r '{\Te c0Sl + dx) d ™+-dI dt j> 



and, if the expression on the right side of this formula be computed, it will be 
found equal to 

{sin (v — V)di — cos (v — P) sin i rfN} X sin N sin i ; 

and, by substituting the values of sin (v — P) and cos (v — P), the same quan- 
tity may be thus written, 

. , „r, , . / »,, , . ?nikT , sin N sin /cos? 

{sin (X — N) d. tan i — cos (K — N) tan i d N} X > 1 8 — ; 

which expression is equal to zero in consequence of what was shown in § 4, 

Wherefore we have, 

and similarly, 

■ ■ (II) 



&iN + |fK = 



> 



rfN 



di 



dN + ^di=0. 



It follows from what has been said that the expressions of dx, dy, dz contain 
dv only, and are independent of the differentials of the five elements, a, e, w, 
N, i, which destroy one another and disappear. And further, if in x, y, z we 
substitute for v, its value in terms of the mean motion and the mean anomaly, 

viz. 

v = I + 2 - O, 
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the expressions of doc, dy, dz will contain dZ, only: for dv contains dZ, only, 
and is independent of the differentials of s, e, s. Thus we have 

_, dx ,„ dx 7 7 dy JV dy 7 7 <?2 7 „ <?2 7 , 

dx = rK dZ j = J - i dt, dy-j-^-^dt, dz = I - K dZ j = j I dt. 

It is in these properties that we recognise the principle of the Variation of 
the arbitrary constants. The finite expressions of x, y, z, being the same in the 
immoveable ellipse described by the sole action of the centripetal force of the 
sun, and in the variable ellipse which represents the motion of a disturbed 
planet, they will verify the equations (A), supposing the arbitrary quantities 

constant, and neglecting the disturbing forces. The velocities -r-., j|, -^ are 

the same whether the arbitrary quantities remain constant or vary ; and thus, 
for a moment of time d t, the motion in the invariable ellipse coincides with 
that of the planet in its real path. But, in the next moment of time, the planet 
will quit the periphery of the ellipse supposed to continue invariable ; because 
the forces in that orbit are different from the forces which urge the planet. In 
the immoveable ellipse the forces in the directions of the coordinates are equal 

to ~d¥ i ~dfi' ~dF' tne arD i trar y quantities being constant ; but, in the case of 

the planet, the like forces are equal to the same differentials augmented by 
the variation of the arbitrary quantities, the additions thus introduced being 

equal to the disturbing forces, \h -j^, \h -j-, p -j~. It is in this manner that an 

elliptic orbit, by the variation of its elements, is capable of representing at 
every moment of time both the velocity of a disturbed planet, and the forces 
by which it is urged. 

And generally, when a dynamical problem admits of an exact solution, the 
arbitrary quantities may be made to vary so as not to alter the velocities 

(J "V (1 11 ft, % ' 

dl> ctt> di' anc * ^ e Editions which the variation of the same quantities makes 

,, . ddx ddii ddz .„ „.,,., 

to the expressions -jj?> ~jf, -jp wul represent new forces introduced in the 

problem. By means of this artifice we may estimate the effect of any disturb- 
ing forces, more especially of such as bear an inconsiderable proportion to the 
principal forces, in altering the original motion of the body. This is the prin- 
mdcccxxxii. 2 F 
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ciple of the Variation of the arbitrary constants, a method which has been much 
discussed, and which is now probably exhausted. It originated in the first 
researches on physical astronomy, and has been matured in passing through 
the hands of Euler, Lagrange, Laplace, and Poisson. The labours of these 
great geometers have raised up a general analytical theory applicable to every 
case, and requiring no more than the substitution of the particular forces under 
consideration. Invaluable as are such extensive views, the application of for- 
mulas constructed on considerations of so general a nature, may not always 
be very ready or very direct, and may require much subordinate calculation. 
In important problems it may be advantageous to separate the principles of the 
method from the analytical processes with which they are conjoined, and to 
deduce the solution directly from the principles themselves by attending 
closely to the peculiar nature of the case. 

Distinguishing the two planets by their masses m and m', the symbol R 
stands for a function of x, y, z, the coordinates of the disturbed planet m, and 
of x', y, z', those of the disturbing planet m'. The expressions of these latter 
coordinates will be obtained by marking all the quantities in the values of 
x, y, z, with an accent, understanding that the accented quantities denote the 
same things relatively to the orbit of m', that the unaccented quantities repre- 
sent in the orbit of m. The function R may be transformed in two ways, 
according as we substitute, for the coordinates, one set of values or another. It 
will be changed into a function of the four independent quantities r, v, N, i, and 
of the like four accented quantities of the planet m', by substituting the values 
of the coordinates obtained in the beginning of this section : and in this case, 
for greater precision, the partial differentials of R relatively to r and v will be 

written with parentheses, thus, (-jr) and (-^ ). When the values of x, y, z, 

in terms of the mean motion t, and of the six elements, a, s, e, vs, N, i, and the 
like values of the other coordinates are substituted, R will be a function of 
the mean motions Z, and £', and of the respective elements of the two orbits. 
In this latter transformation, the partial differentials of R will be written, 
as usual, without parentheses. It may not be improper to set down here the 
expressions of such of these partial differentials as we shall have occasion to 
refer to, 



dR 
da 
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\dr ) ' da, \dr ) a' 

dR _ fdR\ dr dv /dR\ dv , p 

ds —\dr) ' dv' de + \dv) ' de> W 

dR _ /dR\ (dr.dr dv\ SdR\ dv 
de \dr ) ' \de~*~ dv' dej ' \dv) ' de' 

dR _ /dR\ /dr dr dv\ /d_R\ dv_ 
dts \dr) ' \dw'dv' dva) * \dv ) ' dm* 



in which expressions, it need hardly be observed, that -r-, -7-, 2 — , refer to 
this value of v , 

V — £ + £ — 0. 

Proceeding now to reduce the differentials of the elements of the variable 
orbit to the forms best adapted for use, we have this formula for the mean 
distance a, 

— ■ = — — 2 fd l R : consequently, — = 2 d' R. 

Now, when x, y, z are transformed into expressions of £and the elements 
of the orbit, it has been proved that dx, dy, dz contain d g only, and are inde- 
pendent of the differentials of the elements : wherefore, the value d'R will be 
found by differentiating R, making £ the only variable, that is, we shall have, 

,,-,., dR , , dR , , dR , dR T „ 

dR =i^ dx +-dj d y+d^ dz=s -ds d & 

But substituting this value, 



a = 



dR > • • (12) 

a + 2fa? < ^dt > ,\ 

The mean motion £ is denned by this equation, d £ = — — • But, we have, 

a T 

2f2 
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Let n 2 = -^i, n being the constant of the mean motion in the primitive ellipse, 
0: then 



a 

when £ = 0: then 



r • • • (is) 



Taking next the semi-parameter h 2 , we have, by equation (3), 
hdh = r 2 (d'R-^dr): 

but d'R = \jfj) dv + \-jf) d r i wherefore, 

hdh=(^).r 2 dv = * 2 JT=7 2 .( d £)dZ. 
In order to find the value of \-j^J, let the expressions of -^j and^ in the 

d i) d w 

formulas (C), be added : then, since it has been shown that -j- g + ^ = 1, we 
get, 

de + ds7~\dr) W + rf»/ + V^f/ : 

and, because r is a function of v — ot j ^ + ^ = ; wherefore, 
^R ^R _ /iR\ 

ds ' d B7 \dv J' 

Further, because s always accompanies £, or which is the same thing, because 

R is a function of £ + s, we have -^y = -^r : 
consequently, 

dR,dR_ /dR\ 

By substituting this value, 

i ^ 

*» = a (1-/2) + 2f#JT=?(j% + ?S) rf & j 



q - y 



(14) 
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the semi-parameter of the primitive ellipse being equal to a(l — e' 2 ), and its 
eccentricity to d. 
The eccentricity is determined by this formula, 

e 2 __ x _ ± . 
a 

by differentiating, 

7 hdh . h 2 da hdh , ., „. dK ,„ 

*<?* = — jT+a- •-^-=--ir + a ( 1 - e2 )¥?^ : 

and by substituting the value of hdh, 

For the variation of the perihelion we have the formula (7), which may be 
written in this manner, 

2 

— hdh = cos (v — *?) rfe -f- esin (v — ttr) «/»: 

and by multiplying all the terms by e, 

e sin(t> — w) . edvs = -—hdh — cos (v — zs)ede: 
and because ede= — + h 2 d'R, 

e sin (v — «r) . e ^ot = ( — -| j hdh — h 2 cos (t> — ot) <Z' R. 

Ftatta,,«o^+(S)ir=^+^.£. (£).«.: 

and, by substituting this value, 

• / \ j /2 e , cos (» ■— ■«•) A a cos (u — ot)\ , , , 
esm(w— ot) .edz* = ^— + -5 ^J . Atf/t 

-T*-Tv C0S ( v - m )Yd7) r dv - 
Now -g . j- = e sin (v — ts) ; and it will be found that the coefficient of hdh is 

equal to, 

(2 + e cos (w — «•)) . e sin 2 (u — w) j ^ w 

aTT^l 5 ) = T • Te • ^in(y- W ) : 
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wherefore, by substituting and dividing all the terms by e sin {v ~ w) 
edw = — . -j- e . hdh — cos (v — ii)\-j~ Y . r 2 dv. 

But hdk = irj^J f 2 dv, and, 

, .1 (dr , dr dv\ 

— cos (v — vr) = — l ->- + ;r . ~r i ■ 
v ' a \de ' dv de/ 

wherefore, 

, (dv(dW\ . (dr dr dv\ (dU\ "1 r*dv 

ed& =\Te\Jv~) + \Te+dv' d~e)\Tr')]-~V"- 

and, because r 2 dvz= a 2 ^/! — e 2 .d%, 



at? = — - .-j— .a^. ............ (16) 

The variation of s, the longitude of the epoch, must be deduced from the 
equation (9), viz. 

a 2 v' l — e 2 / t t \ dx > j j 

— 2— ^ (ds — drs) = — j^de — dw. 

From this rfe may be eliminated by means of the equation (7), viz. 

2 

— hdh = cos (y — w) Je + e sin (v — w) rfsr ; 
and the result will be 



V l -^ e 2 . cos (w — w) /■» * \ _ 
. _ — \ '(de — dvr) = 

h d h — < cos (v — ts) — e sin (v — m) -p 1 oV. 
Now the coefficient of a" to- is equal to 
(1 - e 2 ) ^ cos (»— *) - ~ ; 

wherefore, by multiplying by r 2 ,we get 

a 2 *Jl — e 2 . cos (v — w) (di — dm) — —Ir.^.hdh 

— a{a(l — e 2 ) cos {v — w) — 2re}dta: 



_2 dt) 
»• ' de 
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substitute now the value of d-vs in equation (16), and that of hdh viz. 
a • s /l=*(™)dZ, then 

cos(v-a) {dz-dxs) = -2r£- (-—) d% 

_ f a(l -e 2 )cosp- CT 2j- l^^r. 

But, as appears from the formulas (C), 

dv SdR\ dR , . /dR\ 

Te' {-d^f^-dV+^^^-^X-d^)'- 

wherefore, by substituting and dividing all the terms by cos (v — nr), 

j j a(l — e 2 ) dR,„ (dR\ J9 

and by substituting the value of da, and observing that f -j^) = -j- • — , we 
obtain 

rft = a N A^(^^)S^-2a 2 if.^ (ir) 

If the formulas (C) be multiplied, each by its own differential, and the re- 
spective results be added, it will be found that the coefficients of \-y-) and 
(^— )are each equal to zero, on account of the equations (8) and (10): so 

that we have, 

dR 7 , dR , , dR . , dR , 
- d ^da + ^de+ 1[7 ds + j^drz = 0: 

and this equation will serve to verify the values of da, de, ds, dzr, which 
have been separately investigated. 

It remains to examine whether the values of di and e?N already found 
(equation (6)), can be expressed similarly to the other elements. The three 
quantities N, P, i, or rather the two N and i, since P varies with N, are inde- 
pendent of r and v, and consequently of •£, a, e, s, m : wherefore, by differen- 
tiating the expressions of x, y, % relatively to i, we shall get 

-jj = r sin (v — P) sin N sin i = z sin N, 
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-^ = — r sin (v — P) cos N sin i = — z cos N, 

dz ™ . rsin(A — N) 
^..= rs\n(v — P)cosi = — ,\ . ; . 

Let these expressions be multiplied respectively by -j^, -jr, -j^, and then 
added ; the result will be 

^ = 4^sinN-^cosN} + ^.^i^: 

di \dx dy J dz */ L + & 

and by substituting the values contained in the formulas (B), 

dR /, ■ „\ dR . .. » T . «R ,, -n. TN 

Ti = C 1 + * 2 ) -17 sm <V~N)- ^x^cos (X-N) : 

and, because * cos (X — N) = ^ sin (X — N), 

tfR . rt XT . f„ . 2 .dU dR ds\ 
■j-. = sm (X - N) . | ( 1 + * 2 ) ^ - ^ • ^ j . 

If the equations (11) be multiplied respectively by ^j, -j-, -j^, and then 
added, this result will be obtained, 
dR 7 . . dR 7 - T .. 

dR 
By combining this equation and the value of -jj with the formulas (6), we get, 

7 xt a J <m JV 1 



di = - v r=7* ' slHl ' £5 • rf & J 



1 <R « (18) 



The differentials of the several elements of the orbit of the disturbed planet 
have now been made to depend upon the function R and its differentials rela- 
tively to the elements themselves and to the mean motion £. Upon the cal- 
culations which this transformation requires, which have long ago been car- 
ried as far as human perseverance can well be supposed to go, we do not here 
enter. The variations of the elements of a disturbed planet, in the most per- 
fect form in which they have been exhibited in the latter part of this paper, are 
the result of the repeated labours of Lagrange and Laplace, who, at different 
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times and by different methods, at last succeeded in overcoming the diffi- 
culties of this great problem. 

In this paper the utmost rigour of investigation has been strictly preserved. 
No admission or supposition has any where been made for the sake of simpli- 
fying calculation or of obtaining a result more readily. The procedure that 
has been followed likewise makes it easy to change the form of the differentials 
of the elements of the orbit, as occasion may require. Thus it is obvious from 
the formulas (C), and from other formulas, that the variations of all the ele- 
ments may be expressed by means of the three functions -^p -j^> -j- ; or, 

J T> J T> 7 TJ 

by means of the three -r— > -j-> -j- ; or, by any two of the differentials of R 

relatively to a, e, e, zs, and one of the two, relatively to i and N ; which re- 
mark is useful in the theory of the comets. 

There is this advantage in expressing the differentials of the elements by 
means of the function R, that inspection alone discovers the nature of the 
terms that enter into every formula. But it is not enough to know the form 
of the terms, we must likewise attend to their convergency. In the present 
state of the heavens there is no difficulty in this respect, because the eccentri- 
cities of the planetary orbits, and their inclinations to the ecliptic are found to 
be small, and it is upon the smallness of these quantities that the convergency 
of the series into which R is developed, mainly depends. In the present cir- 
cumstances of the planetary system, the formulas afford the utmost possible 
facility for computing the inequalities of the elliptic elements. After all, the 
inquiry is difficult enough when it is carried beyond a first approximation ; for 
in the second stage of the process every element that enters into a formula 
being itself a collection of sines or cosines, it is not easy to be assured of the 
nature of the quantities arising from the combination of so many complex 
expressions. 

If we extend our views and consider the stability of the system of the world, 
it is necessary to begin with establishing the convergency of the terms into 
which R is expanded. The mathematical form of these terms will always be 
the same ; but unless their total amount can be estimated with sufficient ex- 
actness by a limited number of them, the human understanding can come to 
no solid decision. Now this will depend upon the effect of the perturbations 
in changing the eccentricities and the inclinations of the orbits to the ecliptic. 

mdcccxxxii. 2 G 
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If it can be proved that these elements, after an indefinite lapse of time, will 
remain of inconsiderable magnitude as they are at present, the convergency of 
the series will be established, and the form of the terms of which R consists, 
will enable us to compute the changes in all the elliptic elements, and to de- 
cide the great question of the stability of the system. But we cannot enter 
upon any extended discussion of these points, and shall conclude this paper 
with some remarks in illustration of the problem we have solved, and of the 
manner in which we have solved it. 

6. If we suppose that there is no disturbing force, or that R = 0, we shall 
have by the equation (I), 

J_ __ JL dx * ± d y* + dz * 

a r ~ ' df 1 . jt* ' 

and if V represent the velocity of the planet at the extremity of r, then, 

V 2 = -r£- , and — = V 2 . 

dr.(t ' a r 

This last equation shows that the mean distance a of an elliptic orbit depends 
only upon the radius vector drawn to any point, and upon the velocity at that 
point. Conceive that the straight line r extends from the sun in a given 
direction and to a given length, and from its extremity suppose that a planet 
is launched into space with the velocity V, the foregoing equation will deter- 
mine the mean distance a of the immoveable ellipse in which the planet will 
revolve. The point from which the planet is projected, and consequently r, 

remaining the same, — and V 2 will vary together ; and if we suppose that a 

becomes equal to a, at the same time that V 2 is changed into V 2 + d . V 2 by 
forces which act continually but insensibly, we shall have these equations, 

d.— = -d.V% and - = --fd. V 2 . 

a ' a a J 

It has been shown that the disturbing forces acting in the directions of x, y, z, 

and tending to increase these lines, are respectively -j^, -j~, -^: and, by the 

principles of dynamics, double the sum of the products of these forces, each 
being multiplied by the element of its direction, is equal to the change effected 
on the square of the velocity : wherefore, 
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and consequently, 

d.±- = -2d'R, andv = 4--2/VR- 

These results agree with the investigation in the fourth section of this paper ; 
and they coincide with the remarkable equation first discovered by Lagrange, 
from which he inferred the invariability of the mean distances and the periodic 
times of the planets, when the approximation is extended to the first power 
only of the disturbing force. _____ 

It has already been observed that -5— . ^— + s is the disturbing force per- 
pendicular to the plane passing through the sun and the coordinate z, that is, 

d R 1 + s 9 
to the planet's circle of latitude ; and likewise that -7— • — — is the disturb- 
ing force in the same plane perpendicular to r the radius vector. The elements 
of the direction of these forces are respectively - — = and y+& : wherefore, 
^R ,. . dR 
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is the variation produced in the square of the velocity in the direction perpen- 
dicular to r. But dv being the small angle described round the sun in the 
time d t, the space described by the planet perpendicular to r, is r dv; and 

T d V 

consequently - is the planet's velocity in that direction. Wherefore, 

CL T> \r M» 

using the symbol i to denote a variation caused by the disturbing forces per- 
pendicular to the radius vector, and observing that these forces produce no 
momentary increase or decrease of that line, we get, 

dv* 



(dR dR \ r*dv* „v 

KdZ dK + -17 ds ) = * • dWTv. = r l 



df .p. 



consequently, 

„ o (dR , , dR , \ , v dv* v. r*dv 9 

2r2 \dT dx + -di ds ) =rAl •df7^ ==l -JWT ( , : 

7^ dv 

and, as - and its square vary by no other cause but the action of the 

CL V \ W> 

forces perpendicular to r, we have 

2r2 (ll dx + d dJ ds ) = «*•_££• 
Now this is the same differential equation that has already been obtained by 

2 62 
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a different method in equation (3) of the second section, and from which the 
value of h 2 , the semi-parameter of the variable elliptic orbit, was deduced. 
That element is therefore as much an immediate deduction from the disturbing 
forces, as is the mean distance in the equation of Lagrange. As the variation 
of a is the effect of the disturbing force in altering the velocity in the orbit, so 
the variation of h 2 is the effect of that part of the disturbing force which alters 
the exact proportionality to the times of the areas described round the sun. 
The two elements are together sufficient for determining both the form and 
the magnitude of the momentary elliptic orbit. The placing of this ellipse so 
as to be in intimate contact with the real path of the planet, a procedure 
which corresponds to finding the relation between the arcs 6 and v, determines 
the motion of the line of the apsides. 

If, lastly, we attend to that part of the disturbing force which is perpen- 
dicular to the circle of latitude passing through the planet, and proceed as 
before, we shall obtain the differential of the equation (4) in the second section. 
This differential is therefore the effect of the disturbing force in altering the 
momentary area which is described in the immoveable plane of xy, and which, 
without the action of this force would be proportional to the time. The elemen- 
tary area in the immoveable plane is the projection of the area described in the 
same time in the plane of the orbit ; the proportion of the two determines the 
cosine of the inclination of the variable plane in which the planet moves ; and 
from this it is easy to determine the position of the line of the nodes, as has 
been fully explained. 

What has been said is independent of the nature of the forces in action ; and 
it is obvious that the same method may be applied to estimate the effect of any 
extraneous force in disturbing the elliptic motion of a planet. 

It would appear that in the view we have taken of this problem, we have 
been making an approach to some general hints contained in the corollaries of 
the seventeenth proposition of the first book of the Principia. A connexion 
between the most recondite results of modern analytical science, and the 
original ideas thrown out by an author who, although he accomplished so 
much, has unavoidably left much to be supplied by his successors, is un- 
doubtedly worthy of being remarked, and may suggest useful reflections. 

December 22, 1831. 



